0.1. Failure Theories 

In the previous section, we introduced the concept of stress, strain and the relationship 
between stresses and strains. We also discussed failure of materials under uniaxial state 
of stress. Failure of engineering materials can be broadly classified into ductile and brittle 
failure. Most metals are ductile and fail due to yielding. Hence, the yield strength 
characterizes their failure. Ceramics and some polymers are brittle and rupture or fracture 
when the stress exceeds certain maximum value. Their stress-strain behavior is linear up 
to the point of failure and they fail abruptly. 

The stress required to break the atomic bond and separate the atoms is called the 
theoretical strength of the material. It can be shown that the theoretical strength is 
approximately equal to E/3 where, E is Young’s modulus. 1 However, most materials fail 
at a stress about one-hundredth or even one-thousandth of the theoretical strength. For 
example, the theoretical strength of aluminum is about 22 GPa. However, the yield 
strength of aluminum is in the order of 100 MPa, which is l/220 ,h of the theoretical 
strength. This enormous discrepancy could be explained as follows. 

In ductile material yielding occurs not due to separation of atoms but due to sliding of 
atoms (movement of dislocations) as depicted in Figure 1.1. Thus, the stress or energy 
required for yielding is much less than that required for separating the atomic planes. 
Hence, in a ductile material the maximum shear stress causes yielding of the material. 

In brittle materials, the failure or rupture still occurs due to separation of atomic 
planes. However, the high value of stress required is provided locally by stress 
concentration caused by small pre-existing cracks or flaws in the material. The stress 
concentration factors can be in the order of 100 to 1,000. That is, the applied stress is 
amplified by enormous amount due to the presence of cracks and it is sufficient to 
separate the atoms. When this process becomes unstable, the material separates over a 
large area causing brittle failure of the material. 



Figure 1.1: Material failure due to relative sliding of atomic planes 


Although research is underway not only to explain but also quantify the strength of 
materials in terms of its atomic structure and properties, it is still not practical to design 
machines and structures based on such atomistic models. Hence, we resort to 
phenomenological failure theories, which are based on observations and testing over a 
period of time. The puipose of failure theories is to extend the strength values obtained 


1 T.L. Anderson, Fracture Mechanics - Fundamentals and Applications, Third Edition, CRC 
Press, Boca Raton, FL, 2006. 
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from uniaxial tests to multi-axial states of stress that exists in practical structures. It is not 
practical to test a material under all possible combinations of stress states. In the 
following, we describe some well-established phenomenological failure theories for both 
ductile and brittle materials. 

Strain Energy 

When a force is applied to a solid, it deforms. Then, we can say that work is done on the 
solid, which is proportional to the force and deformation. The work done by applied force 
is stored in the solid as potential energy, which is called the strain energy. The strain 
energy in the solid may not be distributed uniformly through out the solid. We introduce 
the concept of strain energy density, which is strain energy per unit volume, and we 
denote it by U 0 . Then the strain energy in the body can be obtained by integration as 
follows: 


U = Iff U o( X ’V’ Z ) dV 

V 


( 0 . 1 ) 


where the integration is performed over the volume V of the solid. In the case of uniaxial 
stress state stain energy density is equal to the area under the stress-strain curve (see 
Figure 1.2). Thus, it can be written as 


U 0 = ^ ae (0.2) 

For the general 3-D case the stain energy density is expressed as 

TT 1 

(^0 2 “ 1 “ ^z^z "t" x yz1yz T~zxlzx ^xy^fxy (0-3) 

If the material is elastic, then the strain energy can be completely recovered by unloading 
the body. 

The strain energy density in Eq. (0.3) can be further simplified. Consider a coordinate 
system that is parallel to the principal stress directions. In this coordinate system, no 
shear components exist. Extending Eq. (0.3) to this stress states yields 

U o ~ 2 a 3 £ 3' ^ 


(0.4) 



Chap. 2 Uniaxial Bar and Truss Elements 


5 



From Section 1.3, we know that stresses and strains arc related through the linear elastic 
relations. For example, in case of principal stresses and strains, 

% = ^(^1 i ' a 2 ~ 

e 2 = ^(a 2 - vcr\ - vcr-i) (0.5) 

£ 3 = ^(03 - VOl ~ 


Substituting from Eq. (0.5) into Eq. (0.4), we can write the strain energy density in terms 
of principal stresses as 


U n = 


2 E 


+ of + ^3 “ + a 2 a 3 + a i a 3) 


( 0 . 6 ) 


The strain energy density can be thought of as consisting of two components: one due 
to dilation or change in volume and the other due to distortion or change in shape. The 
former is called dilatational strain energy and the latter distortional energy. Many 
experiments have shown that ductile materials can be hydrostatically stressed to levels 
beyond their ultimate strength in compression without failure. This is because the 
hydrostatic state of stress reduces the volume of the specimen without changing its shape. 


Decomposition of Strain Energy 

The strain energy density at a point in a solid can be divided into two parts: dilatational 
strain energy density Uu that is due to change in volume, and distortional strain energy 
density, U d , that is responsible for change in shape. In order to compute these 
components, we divide the stress matrix also into similar components, dilatational stress 
matrix, a h , and deviatoric stress matrix, o d . For convenience we will refer the stresses to 
the principal stress coordinates. Then the aforementioned stress components can be 
derived as follows: 
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The dilatational component o h is defined as 

_ (Ji ~\~ <J2 + O3 _ O'xx + a yy + a zz 

° h ~ 3 “ 3 


(0.7) 


( 0 . 8 ) 


which is also called the volumetric stress. Note that 3o h is nothing but the first invariant I\ 
of the stress matrix in Eq.. Thus, it is independent of coordinate system. One can note that 
ct/, is a state of hydrostatic stress and hence the subscript /? is used to denote the 
dilatational stress component as well as dilatational energy density 

The dilatational energy density can be obtained by substituting the stress components 
of the hydrostatic stress state in Eq. (0.8) into the expression for strain energy density in 
Eq. (0.6), 


U h =^ol+al+al -2 v{o h a h + a h a h + a h a h ) 


3 (1 - 2 u) 2 
2 E ° h 


and using the relation in Eq. (0.8), 

3 (1 - 2v) eq + <t 2 + 0-3 2 


(0.9) 


U h = 


E 3 

o \ + 0 2 + 03 + 2(<7i<7 2 + 02°3 T 0103) 


1 2 V 2 , 2 1 „2 

6 E " 


( 0 . 10 ) 


Distortion Energy 

The distortion paid of the strain energy is now found by subtracting Eq. (0.10) from Eq. 
( 0 . 6 ), as 


U d =u 0 - u h 


1 + V 2 1 2 1 2 

a l + 02 + 03 — 0102 — 0203 — 0103 

1 + v 01 — 0 2 2 + 02 03 ^ T 03 — 01 

~3E 2 


( 0 . 11 ) 


It is customary to write U d in terms of an equivalent stress called von Mises stress a VM as 

( 0 - 12 ) 


rr _ 1 + v 2 
u d ~ om °VM 


3 E 


The von Mises stress is defined in terms of principal stresses as 
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a VM 


(7i — <T 2 + (7o — CT 3 + CT3 — <7i 


(0.13) 


Distortion Energy Theory (Von Mises) 

According to the von Mises’s theory, a ductile solid will yield when the distortion energy 
density reaches a critical value for that material. Since this should be true for uniaxial 
stress state also, the critical value of the distortional energy can be estimated from the 
uniaxial test. At the instance of yielding in a uniaxial tensile test, the state of stress in 
terms of principal stress is given by: o\ - o Y (yield stress) and er 2 = 03 = 0. The distortion 
energy density associated with yielding is 

Ud = 1 -J^ a Y (°‘ 14 ) 


Thus, the energy density given in Eq. (0.14) is the critical value of the distortional energy 
density for the material. Then according to von Mises’s failure criterion, the material 
under multi-axial loading will yield when the distortional energy is equal to or greater 
than the critical value for the material: 


1 + iy 
3 E 


2 

a VM 


> 


1 + ;y 2 
3 E Y 


a VM > a Y 


(0.15) 


Thus, the distortion energy theory can be stated that material yields when the von Mises 
stress exceeds the yield stress obtained in a uniaxial tensile test. 

The von Mises stress in Eq. (0.11) can be rewritten in terms of stress components as 




{<7xx ~ <7yy) + ( <7yy ~ °zzf + ( a zz ~ °xxf + 6(ri + tL + T^) 


(0.16) 


For a two-dimensional plane stress state, <73 = 0, the von Mises stress can be defined in 
terms of principal stresses as 


a VM 


l 2 2 

V< 7 1 — <J\<i2 + cr 2 


(0.17) 


and in terms of general stress components as 


&VM 




+ <7, 


yy 


^xx^yy A 3T 


xy 


(0.18) 


The two-dimensional distortion energy equation (0.18) describes an ellipse, which when 
plotted on the rT|—rr 2 plane as shown in Figure 1.3. The interior of this ellipse defines the 
region of combined biaxial stress where the material is safe against yielding under static 
loading. 
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Consider a situation in which only a shear stress exists, such that o x = a y - 0 , and x xy 
= z. For this stress state, the principal stresses are o\ - -02 - r and <t 3 = 0. On the 01-02 
plane this pure shear state is represented as a straight line through the origin at -45° as 
shown in Figure 1.3. The line intersects the von Mises failure envelope at two points, A 
and B. The magnitude of o\ and 02 at these points can be found from Eq. (0.17) as 

Oy = erf + er^i + of = 3of = 37^ 


Tmax = <?! = = 0.577<7y 


(0.19) 


Thus, in pure shear stress state, the material yields when the shear stress reaches 0.577 of 
a Y . This value will be compared to the maximum shear stress theory described below. 



Figure 1.3: Failure envelope of the distortion energy theory 


The Maximum Shear Stress Theory (Tresca) 

According to the maximum shear stress theory, the material yields when the maximum 
shear stress at a point equals the critical shear stress value for that material. Since this 
should be true for uniaxial stress state, we can use the results from uniaxial tension test to 
determine the maximum allowable shear stress. The stress state in a tensile specimen at 
the point of yielding is given by: o\ - a Y , 02 = 03 = 0. The maximum shear stress is 
calculated as 


Anax 



. cry 
* TY = ~2 


( 0 . 20 ) 


This value of maximum shear stress is also called the yield shear stress of the material 
and is denoted by x Y . Note that x Y - oy/2. Thus, the Tresca’s yield criterion is that yielding 
will occur in a material when the maximum shear stress equals the yield shear strength. 
x Y , of the material. 

The hexagon in Figure 1.4 represents the two-dimensional failure envelope 
according to maximum shear stress theory. The ellipse corresponding to von Mises’s 
theory is also shown in the same figure. The hexagon is inscribed within the ellipse and 
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contacts it at six vertices. Combinations of principal stresses o x and oi that lie within this 
hexagon arc considered safe based on the maximum shear stress theory, and failure is 
considered to occur when the combined stress state reaches the hexagonal boundary. This 
is obviously more conservative failure theory than distortion energy theory as it is 
contained within the latter. In the pure shear stress state, the shear stress at the points C 
and D correspond to 0.5ov, which is smaller than 0.577 a Y according to the distortion 
energy theory. 



Maximum Principal Stress Theory (Rankine) 

According to the maximum principal stress theory, a brittle material ruptures when the 
maximum principal stress in the specimen reaches some limiting value for the material. 
Again, this critical value can be inferred as the tensile strength measured using a uniaxial 
tension test. In practice, this theory is simple, but can only be used for brittle materials. 
Some practitioners have modified this theory for ductile materials as 

<Ti > a v (0.21) 

where o\ is the maximum principal stress and op the ultimate strength described in. 
Figure 1.5 shows the failure envelope based on the maximum principal stress theory. 
Note that the failure envelopes in the first and third quadrants arc coincident with that of 
the maximum shear stress theory and contained within the distortion energy theory. 
However, the envelopes in the second and fourth quadrants are well outside of other two 
theories. Hence, the maximum principal stress theory is not considered suitable for 
ductile materials. However, it can be used to predict failure in brittle materials. 
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Figure 1.5: Failure envelope of the maximum principal stress theory 



